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ABSTRACT 

We show that the UV spectrum (1280-3200 A) of the "superficially normal" A-star Vega, as observed 
by the IUE satellite at a resolution comparable to the star's rotational broadening width, can be fit 
remarkably well by a single-temperature synthetic spectrum based on LTE atmosphere models and 
a newly constructed UV line list. If Vega were a normal, equator-on, slow-rotating star, then its 
spectrum and our analysis would indicate a temperature of T c e ~ 9550 K, surface gravity of log g 
~ 3.7, general surface metallicity of [m/H] ~ —0.5, and a microturbulence velocity of v tur b — 2.0 



kms 



Given its rapid rotation and nearly pole-on orientation, however, these parameters must 



be regarded as representing averages across the observed hemisphere. Modeling the complex UV 
line spectrum has allowed us to determine the specific surface abundances for 17 different chemical 
elements, including CNO, the light metals, and the iron group elements. The resultant abundance 
pattern agrees in general with previous results, although there is considerable scatter in the literature. 
Despite its peculiarities, Vega has turned out to provide a powerful test of the extent of our abilities 
to model the atmospheric properties of the early A-stars, particularly the detailed UV line spectrum. 
The value of the measurements from this pilot study will increase as this analysis is extended to more 
objects in the rich high-dispersion IUE data archive, including both normal and peculiar objects. 

Subject headings: stars: atmospheres — stars:abundances — stars:individual(Vega) 



1. INTRODUCTION 

The main sequence A-stars (including the range from 
late-B through early-F) include a fascinating array of 
distinct, chemically peculiar objects. Among the spe- 
cific types are the Am stars, HgMn stars, He- weak stars, 
and t he A Bootis stars (see, e.g., the discussion by IWolfil 
1983). The lack of strong stellar winds in this spectral 
range, combined with the absence of a vigorous convec- 
tive zone, results in a quiet outer atmosphere in which 
relatively slow diffusive processes can operate, poten- 
tially producing both vertical and horizontal chemical 
inhomogeneities in the line-forming region of the stars' 
photospheres. A number of specific processes have been 
suggested as contributing to the various forms of pecu- 
liarities, including magnetic fields, rotation, gravitational 
settling, accretion from interstellar space, and accretion 
from circumstellar disks. Indeed, a recent IAU sympo- 
sium was entirely devote d to s tudy of the " A- Sta r Puz- 
zle" (jZverko et all 120041 ) . and iLandstreetl (|2004l ) notes 
how studies of A-stars can lead to insights into the phys- 
ical proce sses pro d ucing the observed phen omena. As 
noted by iCowlevI (|1980t also iCowlev et all 119821) . the 
general diffusion model for the chemically peculiar stars 
holds that they evolve from stars with normal chem- 
ical composition. Likewise for the case of accretion. 
The importance of "normal" A-stars is thus that they 
set the compositional benchmarks against which pecu- 
liar abundance patterns are measured. Even the nor- 
mal A-stars, however, pose challenges — particularly 
in identifying exac tly what normal is fo r an A-star. 
Over 40 years ago iConti &; Stroml (|1968l and see ref- 
erences within) detected compositional anomalies in A- 
stars classified as normal. These stars have since be- 
come known as the "superficially normal" stars, i.e., 



stars whose peculiarities are not recognized until quan- 
titative spectral analyses are performed. The poster 
child for such stars is Vega (HD 172167), the photo- 
metric and spectral classification standard which, nev- 
ertheless, is now well known to be a pole-on, rapidly 
rotating star with extreme gravity darkening, decid- 
edly non-solar surface ab undances, an infrared excess , 
and a dusty disk ( e .g.. iSadakane fc Nishimural 1 1981 ; 
l Aumann et all Il984t iGravl 119851: [ Gulliver et all 11994 ; 
iWilner et all l2002t iPeterson et all 120061 ). Such results 
lead to the question of whether tr uly "normal" A-stars 
even exist ([Hill fc Landstreetlfl99l ! 

This paper represents a pilot study to determine the 
feasibility of modeling the full UV spectrum of the A- 
stars (both normal and peculiar) in order to extract in- 
formation on the basic stellar properties (e.g., T c g and 
log g) and the atmospheric surface abundances for a wide 
variety of chemical elements. The ultraviolet, short- 
ward of ^3200 A, is an extremely rich spectral region 
for such a study, as it contains thousands of individual 
absorption lines from many chemical elements in mul- 
tiple stages of ionization, generally including the dom- 
inant stages. Of course, numerous studies of A-star 
abundances have been performed in the past including, 
for example, the wor k by Adelman and c ollab o rators 
(|Adelman et all I2006L a nd earlier paper s), iHill (|1995l 
and earlier papers), and iVarennd (|1999l ). These stud- 
ies generally involve short segments of high-resolution 
ground-based spectra, one to ~20 elements, and dozens 
to (occasionally) a few hundred individual spectral fea- 
tures. Only a ver y few quantit ative studies have utilized 
UV spectr a (e.g... ISmithlll997l a nd earlier papers in the 
series; and lAdelman et al.ll2004l) and these have concen- 
trated on small segments of data containing spectral fea- 
tures for a small number of elements. The new aspect 
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of this program is that we plan to utilize essentially all 
the spectral information present in the wavelength range 
^1200-3200 A. This study will complement the previous 
work, adding a valuable check on results derived from 
the optical spectral range, adding access to additional el- 
ements, and potentially revealing heretofore undetected 
behaviors in the surface composition of the A-stars. 

The chief data set for this study is to be the archive 
of high-dispersion observations made b y the Interna- 
tional Ultraviolet Explorer {IUE) satellite (Boggess ct al. 
[l978 b.aV These data cover the wavelength range 1150- 
3200 A with a spectral resolution of ~20 kms -1 , utiliz- 
ing observations with a short wavelength camera (1150- 
2000 A) and a long-wavelength camera (1900-3200 A). 
A search of the IUE database, managed by the Multi- 
mission Archive at STScI (MAST), reveals -300 stars in 
object classes 22 (B6-B9.5 V-IV) and 30 (A0-A3 V-IV) 
for which data from both the short- and long-wavelength 
IUE cameras are available. Of these, ~100 are nominally 
normal stars, and the rest exhibit some compositional 
anomaly. The great majority of these spectra have never 
been fully analyzed and thus this archive constitutes a 
treasure trove of virtually untapped information on the 
surface properties of the A stars. 

The specific goals of this planned study are: 

1. Determine the degree to which the full UV spec- 
trum can be reproduced using the standard tools 
of LTE atmospheric models and LTE spectral syn- 
thesis; 

2. Search for systematic discrepancies which may sig- 
nal shortcomings in the input physics or assump- 
tions, or both, of the standard tools or, alterna- 
tively, specific peculiarities with the target stars; 

3. Examine the completeness and accuracy of the line 
lists used to generate the synthetic spectra; 

4. Develop better estimates of physical properties of 
the target stars and, particularly, confirm, refine, 
and expand surface abundance determinations us- 
ing a different technique and a generally different 
set of spectral features than most previous analy- 
ses; 

5. Examine the abundance patterns revealed from the 
study of a large group of normal and peculiar A- 
stars to gain insight into the processes driving the 
compositional anomalies. 

The target for this pilot study is the aforementioned 
Vega. Given its special nature, its use for the pilot 
study might be questi oned. Echoing the thoughts of 
iGarria-Gil et al.l ()2005l ). however, we note that the qual- 
ity and quantity of data available for Vega are unique. 
We will be able to compare our results to a large number 
of previous studies — based on similar physical assump- 
tions — to help determine the reliability of the current 
analysis, while perhaps adding additional insights into 
this interesting star. This paper is not expected nor in- 
tended to supersede other studies which more fully take 
into account Vega's particular physical peculiarities, no- 
tably, the high rotation rate which is partially m asked 
by a near pole-on orientation (e.g.. iHill et al.ll2010D . In 



concert with such studies, however, this analysis will help 
in determining the limitations in accuracy which may be 
expected when simplified, though still computationally- 
challenging, models are used in the analysis of complex 
objects. 

The outline of this paper is straightforward. In fj2]we 
describe the data sets assembled for this study, including 
the IUE data noted above, as well as a limited number 
of UV spectra from other instruments. The modeling 
process, through which we attempt to reproduce the full 
UV spectrum, is presented |J3J The results of the mod- 
eling process are described in 2] and §5\ giving both a 
closeup of the detailed spectrum and a broad view of 
the spectral energy distribution, respectively. The sur- 
face composition of Vega's atmosphere implied by our 
analysis is discussed in fJBJ along with a comparison with 
previous results. Some final comments are given in [J7J 

2. THE DATA 

The primary data set used here consists of IUE spec- 
trophotometry obtained using the large-aperture, high- 
resolution echelle observing mode with both the short- 
wavelength SWP camera (1150-2000 A) and the long- 
wavelength LWP camera (1900-3200 A). The spectral 
resolutions of the SWP and LWP data, ~23 kms" 1 and 
~16 kms -1 , respectively, are comparable to the rota- 
tional broadening observed for Vega (see ©. The spe- 
cific IUE spectra used are listed in Table [T] The six 
SWP observations were all obtained during the time pe- 
riod 1988-1992 and have similar exposure times and data 
quality. Two other large-aperture, high-resolution SWP 
spectra exist, from earlier time periods and with differ- 
ent exposure times, but were excluded in the interests of 
uniform data quality. The five LWP observations listed 
in Tabled] were obtained between 1986 and 1992 and all 
have identical exposure times. One other LWP spectrum, 
from an earlier time period and with a shorter exposure 
time, was excluded for the sake of data uniformity. For 
the same reason, three high-resolution spectra from the 
LWR camera, which has a different photometric response 
and a different spectral resolution from the LWP camera, 
were not included. 

The 11 IUE spectra in Table [T] were all acquired from 
the IUE data archive as hosted by MAST0. These spectra 
were processed with the New Spectral Image Processing 
System (NEWSIPS; INichols fc L inskvl Il996l ) . We pro- 
duced a single mean SWP spectrum by averaging the six 
individual SWP spectra and, likewise, produced a sin- 
gle mean LWP spectrum. The averaging was done using 
routines in the IUE Data Analysis Center (IUEDAC) 
software suite, written in the Interactive Data Language 
(IDL) and downloaded from MAST. In particular, we 
used the routine CRSCOR to measure any small wave- 
length shifts that might exist within the sets of SWP and 
LWP spectra and then used COADD to correct the shifts 
and average together the individual spectra. The aver- 
aging was done on an echelle order-by-order basis, with 
the mean SWP spectrum consisting of 51 separate orders 
(covering the wavelength range 1151-1978 A) and the 
LWP spectrum consisting of 47 orders (covering 1899- 
3097 A). The final step in the data processing was to trim 

1 The MAST Archive was accessed at 

http://archive.stsci.edu/hst/ 
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the spectra to eliminate low-quality data and overlap be- 
tween orders. Specifically, we: 1) removed the first and 
last 10% of the data points from each order, 2) removed 
points which overlapped with the next order (in the sense 
of increasing wavelength), 3) eliminated all SWP data 
shortward of 1282 A (where the noise level is high), and 
4) eliminated all LWP data shortward of 1979 A (which 
overlaps with the higher quality SWP data). The re- 
sult of this procedure is nearly complete coverage of the 
UV between 1282 and 3097 A, with a spectral resolution 
comparable to the expected intrinsic width of the stel- 
lar features and consisting of data from 80 different IUE 
echelle orders (38 from the SWP and 42 from the LWP). 
This spectrum is shown in the top panel in Figure[TJ The 
apparent noisiness of the spectrum is mainly the presence 
of many strong absorption features, although the short 
wavelength end of the LWP data (1980-2300 A) does have 
a significantly lower signal-to-noise ratio than other re- 
gions of the spectrum. Small gaps at the long wavelength 
ends of the SWP and LWP data arise from incomplete 
overlap between adjacent echelle orders. Some medium 
frequency structure in the data, most noticeable around 
2200 A, arises from incomplete removal of the echelle 
"ripple," a long-standing problem with high-resolution 
IUE data. The analysis described in $3] below compen- 
sates for this effect. The model spectrum in the bottom 
panel of Figure [1] will be discussed in 

Our analysis also utilizes a small number of high qual- 
ity, medium- and high-resolution UV spectra of Vega ob- 
tained by the Goddard High Resolution Spectrograph 
(GHRS) aboard the Hubble Space Telescope (HST). 
Likewise, we include a far-UV spectrum obtained by 
the Berkeley Extreme and Far-Ultraviolet Spectrome- 
ter (BEFS). The BEFS instrument was flown aboard 
the ORPHEUS telescope on the ORPHEUS-SPAS space 
sh uttle missions in 19 93 and 1996 and is described fully 
bv lDixon etaLl (pool . The relevant GHRS and BEFS 
observations are listed in Table [1] All these data were 
obtained from MAST and no additional processing was 
performed. 

Finally, low-resolution UV spectra from the IUE satel- 
lite provide a reference for the absolute flux levels 
of Vega, which will be examined in |J5] To facili- 
tate this discussion, we obtained from MAST a set 
of 7 low-resolution SWP spectra, consisting of im- 
ages SWP27024, SWP29864, SWP29866, SWP29867, 
SWP30548, SWP30549, SWP32906, and SWP32907. 
These data were all acquired by IUE using the "fast- 
trail" technique developed for bright sources and the 
list includes all the SWP spectra for Vega for which 
no detector saturation occurred. We corrected the 
data — which were processed using the NEWSIPS sys- 
tem — for residual systematic calibration err ors and 
placed them onto the HST/FOS flux scale of iBohlinl 
(1996) using the corr e ctions and algorithms described by 
iMassa fc Fitzpatrickl (|2000L hereafter MF00). The indi- 
vidual spectra were then averaged to form a single far-UV 
spectrum covering the wavelength range ~1250-2000 A 
and with a spectral resolution of approximately 6 A. 

3. THE ANALYSIS 

Our analysis consists of computing a rotationally- and 
instrumentally-broadened LTE synthetic UV spectrum 



from a pre-existing stellar atmosphere model and com- 
paring it with the observed high-resolution UV spectrum 
of Vega. Adjustments are then made iteratively in the 
various stellar properties which determine the synthetic 
spectrum until the best possible agreement is obtained 
with the observations. For t he atmo s phere models, we 
use the ATLAS9 grid from iKuruczl (|1991[ ). computed 
with a microturbulence velocity of 2 kms -1 and a wide 
range of T e e, logg, and [m/H] values, spanning the re- 
gion of interest. To compute the synthetic spectrum, we 
use the program SPE CTRUM (Version 2. 76e ) developed 
by Richard Gray fsee lGrav fc CorballvlfiOQlR . 

The main analysis was performed on the high- 
resolution IUE data described in the previous section, 
which consist of 80 separate echelle orders covering the 
range 1282-3097 A, with no overlap between orders and 
all normalized by the mean flux as computed over the 
entire region. The analysis includes a set of global stel- 
lar parameters and a set of order-specific parameters. 
The former include the effective temperature T c fj, sur- 
face gravity logg, and overall abundance [m/H] — all of 
which specify the ATLAS9 atmospheric structure model 
- as well as a set of individual elemental abundances 
X/H, the microturbulence velocity v tur b, and the rota- 
tional broadening velocity v sin i — all of which are used 
in the running of SPECTRUM and the construction of 
the synthetic spectrum. The order-specific parameters 
(i.e., those determined individually for each echelle or- 
der) include the radial velocity v ra d and a 5-point "con- 
tinuum" used to match the low-frequency shape of the 
individual echelle order data and the synthetic spectrum. 
V ra d is taken as an order-dependent quantity to help 
compensate for order-specific wavelength calibration un- 
certainties. The "continuum" is computed using a cu- 
bic spline interpolation between 5 equally-spaced anchor 
points (for each order) whose values are determined in 
the analysis. The use of this continuum removes the in- 
fluence of low- and medium-frequency phenomena on our 
analysis, including interstellar extinction, shortcomings 
in the absolute calibration of the data, and shortcomings 
in the removal of the echelle "ripple" signature from the 
IUE data. Our results therefore are dependent only on 
the high-frequency information in the data, i.e., the spec- 
tral lines, and insensitive to the general level or shape of 
the spectral energy distribution (SED). In the analysis, 
predetermined values can be adopted for any or all of 
the parameters described above, or best-fit values can be 
determined iteratively. 

Beginning with an initial set of estimates for the full 
set of free parameters, the analysis proceeds as follows: 

1. Construct a stellar atmosphere model for the de- 
sired set of T c ff, log g, and [m/H] values. We 
produce this model by logarithmically interpolat- 
ing the structure variables within the ATLAS9 2 
kms -1 grid using the set of models on the vertices 
of the grid which surround the desired parameter 
set. The result is then outputted to a file in the 
format of the ATLAS9 models and serves as input 
for the second step. 

2. Compute the full resolution synthetic UV spectrum 

2 SPECTRUM was obtained from 

http://wwwl.appstate.edu/dept/physics/spectrum/spectrum.html 
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covering the full wavelength range of the data using 
the program SPECTRUM, whose inputs include 
the atmosphere model, a set of elemental abun- 
dances, a value of Ut ur b, and a line list containing 
the atomic data for the relevant set of absorption 
features. 

3. Rotationally broaden the full spectrum, using the 
desired desired value of i>sini. 

4. Slice the rotationally broadened spectrum into 80 
segments covering the wavelength ranges of the in- 
dividual IUE orders. These segments range from 
12 to 34 A wide. Broaden each segment with the 
appropriate instrumental resolution, velocity-shift 
them according to the order-dependent value of 
v ra d, and the scale them by the interpolated 5-point 
"continuum." 

5. Compare this set of theoretical spectra with the 
observed IUE data and make adjustments to the 
free parameters to improve the agreement. 

6. Go back to step 1 and repeat the process until no 
further improvement in agreement between theory 
and observation can be obtained. The final set of 
free parameters constitutes our "best-fit" model. 

The iterative comparison between theoretical and ob- 
served spectra is handled by the non-linear l east squares 
proced ure MPFIT written in the IDL by Markwardt 
(2009). Weighting factors for the analysis were deter- 
mined from measurements of the signal-to-noise ratio 
across each echelle order. Weights were set to zero for 
cases of detector blemishes and in the cores of strong 
resonance lines which might be expected to have some 
interstellar contamination. Also, the weights were set to 
zero for regions including some C I and Si I lines which 
are found to have suspect atomic parameters (see SJSJ> . 

Rotational and instrumental broadening is performed 
with the programs SMOOTH2 and AVSINI which come 
with Gray's SPECTRUM package. Note that the rota- 
tional broadening step is only an approximation, since 
it assumes that Vega is an equator-on star. However, 
departures of the line shapes from the classical rota- 
tional profile seem import ant only for weak lines of 
non-dominant species (e.g. iGulliver et al.l[l99"l . which 
play only a small role in our analysis. The appropri- 
ate FWHM values of the instrumental broadening for 
high-resolution IUE LWP and SWP data were taken 
from Figures 2.21 and 2.23, respectively, of the IUE 
NEWSIPS Information Manual, as obtained from the 
MAST Archives. For the SWP, a value of 3 pixels was 
adopted for all orders, corresponding to 0.10 A (^23 
kms^ 1 ) at the shortest wavelengths and 0.15 A (^23 
kms -1 ) at the longest wavelengths. For the LWP, the 
resolution is somewhat wavelength-dependent, varying 
between 2.2 and 3.0 pixels. This produces a range of 
0.12 to 0.19 A in wavelength space and 16 to 20 kms^ 1 
in velocity space. For each of the 80 IUE spectral or- 
ders, a single value of the spectral resolution, specified in 
A and computed for the mean wavelength of the order, 
was used in SMOOTH2. The broadening profile is as- 
sumed to be Gaussian. Vacuum wavelengths were used 
for all UV calculations. 



Reproducing a high-resolution spectrum requires a re- 
liable list of absorption lines in the spectral range of in- 
terest and an accurate set of atomic parameters for those 
lines. For this study, and for a parallel one of the bright 
A-star Sirius A (Fitzpatrick & McClain, in preparation; 
hereafter Paper II) we constructed a new list, containing 
the most recent updates to the atomic transition knowl- 
edge base. We formed this list by compiling the data 
for all atomic transitions producing spectral lines be- 
tween 1000 A and 3400 A and with absolute line strength 
information (i.e., loggf), utilizing the Vienna Atomic 
Line Database (VALD; iKupka et ani2000f FI and the Na- 
tional Institute of Standards and Technology Atom ic 
Spectra Database (NIST/ASD; iRalchenko etall IMT^ FI. 
The NIST/ASD values of loggf, which represent experi- 
mental determinations, were adopted in cases where they 
and VALD disagreed. This process yielded a list of over 
253,000 features from elements H through U, many with 
radiative, Stark, and van der Waals broadening half- 
widths available. The H lines (Lyman-a and Lyman-/3) 
were removed from the list since they are treated sepa- 
rately by SPECTRUM. Examination of synthetic spec- 
tra produced by this initial list revealed a number of 
clear shortcomings, generally in the form of predicted 
features where no (or only very weak) observed features 
were present, suggesting severe errors in the loggf val- 
ues. We treated these by either removing the predicted 
lines from the list or by empirically measuring log gf us- 
ing the observed strengths of the lines. The details of 
these and other modifications made to the list are given 
in Paper II, where we also provide a table of the em- 
pirical log gf values, which may be of interest to other 
investigators. Likewise, the final line list will be made 
available from the author. As will be noted further in 
the detailed discussion of the derived abundances (ij6]), 
lines with empirically-derived log gf values were masked 
in the main analysis and do not contribute to the results. 

4. THE BEST-FIT MODEL 

After considerable experimentation with the analysis 
technique, we determined that we could solve for robust 
values of T e g, log g, Vturb, v ra d, and the abundances of 17 
chemical elements, while adopting fixed, predetermined 
values for v sin i and the metallicity of the ATLAS9 struc- 
ture models. The results of this analysis are presented in 
Tabled! 

The stellar properties listed in column 1 are mostly 
self-explanatory although a few require some elabora- 
tion. These are: (1) the quantity [m/H], which refers 
to the general metallicity of the ATLAS9 models used to 
provide the atmospheric structure; (2) the radial veloci- 
ties, D ro ,!(SWP) and D ra d(LWP), which are mean values 
for the SWP and LWP cameras, respectively; and (3) 
the individual elemental abundances, e.g., C/H, which 
are the values determined by the analysis and are given 
in the form X/H = log(N x /N H ) + 12.0. The parame- 
ters of our "best-fit" model for the UV absorption line 
spectrum of Vega, as represented by the high-resolution 
IUE SWP and LWP observations, are listed in column 

3 The VALD was acc essed at 

http : / / vald. astro, univie. ac. at/ ~ vald / php /vald.phpl 

The N1ST Atomic Spectra Database was accessed at 
http:/ /physics. nist.gov/asd3. 
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2 of the Table ( "FIT A"). For this fit, we adopted the 
ATLAS9 model grid with roughly 1/3 solar abundances 
(i.e., [m/H] = — 0.5) and assumed a value of usini 
21.x kms- 1 from lGulliver et al.l (fl99l . (The sensitivity 
of the results to both assumptions will be examined be- 
low.) All other parameters in the analysis, including the 
low frequency mapping of the surface fluxes to the ob- 
served fluxes, were free parameters determined by the % 2 
minimization procedure. We fixed vsini due to concerns 
that the IUE spectral resolution was not high-enough 
nor well-enough determined to allow us to defend a new 
measurement of the rotational broadening. In view of 
the insensitivity of the our results to vsini, the adoption 
of a fixed value for the main analysis is reasonable. For 
the same reason, we also made no attempt to model the 
effects of macroturbule nce, whose presenc e has been sug- 
gested by the work of lYoon et alt (|2008l) and Mill et all 
(201(3). As experience with other stars accumulates, we 
will attempt to refine the measurements of the IUE's 
spectral resolution, possibly allowing more information 
to be extracted from the observed line profiles. The fi- 
nal column of Table [5] lists the approximate number of 
spectral features available for the individual abundance 
determinations. Numerous weak features have not been 
included in this tally and, in some cases, the abundances 
are dominated by a smaller number of prominent features 
(see g6]). 

The temperature we find from the UV line spectrum, 
9547 ± 17 K, is generally consistent with results already 
in the literature, mostly based on the level or shape of the 
SED in the optical-throug h-UV spectral region. See Ta- 
ble 3 of IGarcfa-Gil et al.l 12005) for a summary of such 
results, which generally lie in the range 9600-9700 K. 
A different approach, that of modeling the rotationally 
broadened Fe I and Ti II line s at 4 528 A and 4529 A, 
respectively, bv lGulliver et al.l ()1994D . yielded a temper- 
ature of 9695 ± 25 K. Our resul t is pa rticularly close 
to the those of lCastelli fe Kuruczl (I1994I. i.e., 9550 K as 
based on optical colors). ICiardi et al.l (|2001l i.e., 9553 
K as bas ed on the bolometric flux and angular diame- 
ter), and lHill et al.l (|2010l i.e., 9560 K, which represents 
the mean across a rapidly rotating and highly gravity- 
darkened surface). These latter results and ours may 
represent a real convergence on an apparent temperature 
for Vega in the neighborhood of 9550 K . 

The surface gravity in Table [2] 3.72 ± 0.03, is some- 
what lower than most past results. For example, anal- 
yses utilizing the shape of the UV-through-optical SED 
and those using B aimer line profiles g enerally find results 
in the range logq = 3.90-4.00 (e.g.. iCastelli fe Kuruczl 
19941: IFitzpatrick fe Massalll999l: IGarcfa-Gil et al.l 120051: 
Hill et aUhoiOD . On the other hand, iGulliver et all 
( 1994)'s line profile analysis yielded 3.75±0.02, very sim- 
ilar to our result. In considering these values, it should 
be kept in mind that the rapid rotation of Vega, and the 
subsequent extreme gravity-darkening, make the concept 
of a single surface gravity (and a single effective temper- 
ature) invalid. All results for Vega reflect some sort of 
average across the visible surface of the star and it is not 
obvious that different techniques should necessarily yield 
the same results. We will return to this point later in SJ71 

Many analyses of Vega adopt a microturbulence 
broadening velocity of 2 kms -1 , often based on the 



work of Castelli and Kurucz ( e.g.. ICastelli fe Kuruczl 
119941) . IFitzpatrick fe Massal dl999|) found a simi- 
lar value, 1.9 ± 0.4 kms -1 , from modeling the low- 
resolution UV-through-optical SED using ATLAS9 mod- 
els. Results from elemental abundance analyses gen- 
erally yield com patible results, ranging between 1.5 
and 2 kms~He . g .,[Sadakaneet^I|[l98i IQiu et al.ll200TI; 
ISaffe fe Levatol l200ll Our result here, 2.04 ± 0.02 
kms -1 , is thus consistent with th ese previous studies. 
On th e other hand, ana l yses b y lAdelman fe Gulliverl 
(fl990l) . IHill fe Landstreetl fl993l ). and lHill et al.l (|2010f ) 
find values of v tur b closer to 1.0 kms" 1 . We have no ex- 
planation for the discrepancy among these results but, 
clearly, a complete understanding of Vega's atmosphere 
cannot be claimed until they are resolved. 

The radial velocities noted in Table [2] are mean values 
computed from the results for the individual spectral or- 
ders of the IUE SWP and LWP cameras. The uncertain- 
ties listed are actually the sample standard deviations 
for the 2 cameras. The small values, roughly 10% of the 
camera resolutions, testify to the fidelity of the wave- 
length calibrations within each camera. However, a sys- 
tematic offset of more than 17 kms -1 between the two 
cameras is clear, with the SWP result compatible with 
the expected radial velocity of Vega (see, e.g., the set of 
measurements recorded in the SIMBAD database). Our 
treatment of v ra d as a free parameter for each spectral 
order eliminates such problems as factors in the analysis. 
We defer a discussion of the metallicity of Vega, and the 
individual elemental abundances until S}6] below. 

The uncertainties listed for FIT A in Tableware (with 
the exception of those for v ra d as noted above) 1-a sta- 
tistical uncertainties resulting from the % 2 minimization 
procedure. Thus they underestimate the true uncer- 
tainty since they fail to take into account errors in the 
various assumptions built into the analysis and in the 
fundamental physical validity of the model atmosphere. 
In an attempt to gauge the robustness of the analysis, we 
ran a number of fits to the data with significant changes 
in the input assumptions. These results are also shown in 
Tabled in columns 3-10. The tests involved: 1) chang- 
ing the base metallicity of the ATLAS9 models by ±0.5 
dex from that assumed in the best-fit model (FITS B 
and C); 2) changing the assumed value of vsini by ±3 
kms -1 from that in the best-fit model (FITS D and E); 
forcing log g to values ~ ±0.2 dex different than the best- 
fit result (FITS F and G); and 4) forcing Vturb to values 
~ ±0.5 kms -1 different from that found in the best fit 
(FITS H and I). These tests still don't allow us to place 
strict absolute error bars on our results, but do give a 
strong indication that the analysis is very stable against 
uncertainties in values of the assumed parameters. We 
will discuss these results in more detail in 

An overview of the best-fit model can be gained from 
the lower panel of Figure [T] Shown are the surface fluxes 
of the model, broadened by v sin i and the IUE spectral 
resolution, but not mapped to the observed flux levels. 
The chief features of the model, still visible in the com- 
pressed format of Figure [TJ are the steeply declining flux 
in the far-UV (as the wings of the strong H I Lyman- 
a line intrude in the observed wavelength window), the 
1400 A Lyman-a satellite feature, a C I edge at 1446 
A due to ionizations out of the 2s 2 2p 2 level at 21648.0 
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cm 1 , a weaker Si I edge near 1520 A due to ionizations 
out of the ground state, the strong resonance absorption 
lines of Mg II at 2800 A, and many thousands of other 
absorption lines spread throughout the observed spec- 
trum. The 1446 A and 1520 A features are difficult to 
discern in the observed spectrum in the top panel of Fig- 
ure [T] due to calibration deficiencies with the data. The 
1400 A feature arises from photon absorption during an 
interaction between an H° atom and a nearby H + ion. 
It has been seen in the spectra of DA white dwarfs (e.g., 
Greensteml Il980f) . and was first correctly identified by 
Nelan fe Wegnerl (|1985f l and lKoester et al.1 (|1985| ). 

A more detailed view of the best-fit model and its 
level of agreement with the IUE high-dispersion data is 
shown in Figure 12.11 for the three shortest-wavelength 
SWP echelle orders (numbers 105, 106, and 107 cover- 
ing the total wavelength range 1282.22-1318.95 A). The 
model and the data are indicated by the thick and thin 
curves, respectively. Data points excluded from the anal- 
ysis due to detector artifacts are identified with an "x," 
as at 1292.5 A. The full set of Figures 2.1-2.27, which 
cover the entire wavelength range and all 80 SWP and 
LWP echelle orders included in this analysis, are available 
in the online version of the Journal. An ascii file listing 
all the absorption features that contribute to the model 
spectrum can be obtained from the author. A close look 
at these figures shows that, in general, the quality of 
the agreement between the model and Vega's UV line 
spectrum is very good, without any obvious systematic 
discrepancies. The range of lines fitted by the model 
spans weak features, barely visible above the noise level, 
to very strongly damped features. The discrepancies that 
do occur are usually isolated in wavelength space and are 
generally obviously due to spectral features missing from 
the atomic line list or to features for which the available 
atomic data are likely inaccurate, with both issues occur- 
ring more frequently at wavelengths shortward of ^1500 

A. 

The existence of a small number of //ST/ GHRS spectra 
for Vega (as listed in Table [lj allows a more critical ex- 
amination of the agreement between theory and observa- 
tions than possible for the photometrically-inferior IUE 
data — although only for a limited wavelength range. 
This comparison is shown in Figures |3.1H3.4[ for the four 
G160M and G200M observations, and in FigureH for the 
three high-resolution ECH-B observations. For clarity, 
the medium resolution observations in Figures 13.11 13 41 
have been stretched over two panels each. We computed 
the best-fit model exactly as prescribed above, using the 
parameters listed in column 2 of Table [5J We then broad- 
ened the model using the appropriate instrumental res- 
olutions and determined best-fit values of the velocity 
offset (v ra d) and the 5-point "continuum" (used to map 
the model surface fluxes onto the normalized fluxes) for 
each observation. The GHRS G160M and G200M spec- 
tral resolutions were evaluated from the formulae 

(A/AA) G i6om = 16.82A - 2986.4 and (1) 

(A/ AA) G20 oM = 16.20A - 6940.0 , (2) 

obtained from Glenn Wahlgren (private communication) . 
For each observation, the value of AA was evaluated at 
the central wavelength and adopted for the whole obser- 



vation. These FWHM values were 0.070 A (17 kms" 1 ), 
0.069 A (16 kms" 1 ), 0.069 A (16 kms" 1 ), and 0.080 
A (13 kms -1 ) for the observations in Figures 13.21 
13. 3\ and 13.41 respectively. We adopted a Gaussian with 
FWHM = 3 kms" 1 for the ECH-B resolution, as ob- 
tained from the GHRS Instrument Handbook^ 

Figure 3 reinforces the comments above, based on the 
IUE data, that the quality of the fit is good in gen- 
eral, but better at the longer wavelengths. Note that the 
log gf values of the prominent C I lines near 1289 A and 
1310 A, which all arise from the 2s 2 2p 2 level at 10192.63 
cm -1 , were empirically determined by us from the spec- 
tra of Vega and Sirius A (see Paper II) . The values ob- 
tained from VALD were clearly inaccurate and produced 
features much weaker than observed. These lines (indi- 
cated by the open circles in Figure 3) were masked off 
during the current analysis and so did not influence the 
determination of the C abundance. A close inspection of 
model and observations in Figures 13.21 and 13.31 suggests 
that the model spectrum may have been slightly over- 
broadened relative to the data (see, e.g., the region near 
the strong C I lines at 1280 A). Such an effect is not obvi- 
ous in Figure 133] (at least partly due to the noise level) or 
in Figure 13.41 These impressions are seemingly verified 
by allowing the fitting process to choose optimal values 
for v sin i when fitting the GHRS data. The results of 
are 18.0 kms" 1 for the G160M data of Figures l3~2l and 
EH and 20.2 kms" 1 for the G200M data of Figure[32fl 

This issue of line widths can be examined further 
with the three GHRS/ECH-B spectra shown in Figure 
[4j These data allow the most accurate determination of 
vsini, since their spectral resolution (3 kms" 1 ) is sig- 
nificantly narrower than the rotational broadening. The 
only accommodations to the data were the determina- 
tion of velocity offsets and 5-point "continua" to map the 
flux levels for each spectrum. The standard parameters 
clearly provide very good fits to the observed features 
(mostly due to Fe II). If we allow vsini to be deter- 
mined by the ECH-B data themselves, a value of 20.7 
kms" 1 is found, w hich is in excelle nt agreement with 
the recent results of iHill et al.l ((20101 . i.e., vsini = 20.8 
± 0.2 kms" 1 ). This is reasonably consistent with our 
assumed value of 21.8 kms" 1 for the IUE data analysis, 
and broader than suggested by the GHRS G160M spec- 
tra. It may be that the instrumental broadening adopted 
for the G160M data are slightly overestimated, although 
we cannot rule out that this is an artifact resulting from 
inaccuracies in the atomic line list. This issue has no 
bearing on the results of our analysis, which were car- 
ried out with the IUE data. Further, the value of vsini 
suggested by the ECH-B data is comfortably close to the 
assumed value of 21.8 kms" 1 and Table [2] shows that 
our results are insensitive to the exact value adopted for 
v sin i (due to the fact that the spectral resolution of the 
IUE data themselves are comparable to v sin i) . Note the 
presence of interstellar contamination arising from the lo- 
cal ISM in some of the line cores in Figure 0] (indicated 
by the x's) . A discussion of this and a detail of the Fe II 
line at 2344 A can be found in Lallcmcn t et al.l (|1995l) . 

The very high-resolution of the ECH-B data allows for 

5 The GHRS Handbook was accessed at 

http: / / www. stsci. edu / instrument-news/handbooks / ghrs / GHRS_ 1 . html 
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a potentially detailed study of the UV line profiles, akin 
to the attention received by optical line profiles. Such 
studies can involve corrections to line wavelengths or os- 
cillators strengths, or both, and (particularly important 
in the case of Vega) the effects of its high rotation speed 



particularly in their cores (see, e.g., 


Gulliver et al.lll994t 


Yoon ct al. 2008; Takeda ct al. 2008 


). Unfortunately (or 



mercifully), the amount of such data available is far too 
small to contribute significantly to these efforts and the 
comparison in Figure |4] suggests that only mild adjust- 
ments in line position or strengths would be required to 
achieve excellent agreement. The relatively low spectral- 
resolution of our main dataset — the IUE high-resolution 
spectra — and their lower photometric quality limit their 
usefulness for detailed profile analysis. However, the abil- 
ity of the current model to reproduce the degree of de- 
tail available in the IUE high-resolution data, spanning 
nearly 2000 A and including thousands of spectral fea- 
tures, as well as the UV through optical SED, is impres- 
sive and lends confidence in the derivation of elemental 
abundances from the UV spectral synthesis. 

5. THE SPECTRAL ENERGY DISTRIBUTION 

A more global view of our results is shown in Figure 
[S] where we compare the observed SED of Vega (thin 
curve) with the predictions of our best-fit model (thick 
curve) over the range 1050 to 6700 A (thick line). The 
observed SED of Vega was taken from the CALSPEC 
da tabase, file "alpha_lyr_stis_005.ascii," and is discussed 
by iBohlinl (|2007t n. This SED is a composite, based on 
low-resolution IUE SWP data, HST/STIS G230LB and 
G430L data, and Kurucz atmosphere models. We mod- 
ified it in two ways. First, we replaced the spectrum 
below 1200 A (based on a model) with the BEFS data 
described in SJ2] Second, we replaced the data in the 
range 1250-1978 A (based on IUE and HST/STIS) with 
the mean low-resolution IUE SWP spectrum processed 
as described in Sj2] We did this because we noticed a 
large discrepancy in the flux levels between our processed 
IUE spectrum and the CALSPEC fluxes between 1300 
and 1800 A), in the sense that the CALSPEC fluxes 
are larger, in a smooth wavelength-dependent manner 
and by as much as nearly 9% near 1500 A. Outside the 
1300-1800 A range the fluxes agree very well. We have 
never encountered such a discrepancy when modeling 
low-resolution IUE spectra of early-type stars using data 
processed by the MF00 algorithms and it made us suspi- 
cious of the short wavelength end of the CALSPEC SED. 
Hence, we made the wholesale replacement of the 1250- 
1978 A data. (We note that, at longer wavelengths, IUE 
LWP spectra processed with the MF00 algorithms agree 
well with the CALSPEC SED.) Ultimately this modifi- 
cation of the CALSPEC SED has no effect on the con- 
clusions of this study, but its motivation suggests that 
the far-UV fluxes of Vega have not yet been conclusively 
established. 

The model spectrum shown in Figure [5] was computed 
using SPECTRUM and the parameters listed in column 
2 of Table [3J with the spectral range extended outside 

6 The C ALSPEC databa se was ac cessed at 
http: / / archive, stsci . edu/hst / observatory/edbs /calspec . html 



the fitting region, and then broadened to match the 
instrumental resolution of the datasets comprising the 
observed SED. The resulting surface fluxes were then 
scaled to match the flux of the observed SED at 5556 
A, namely 3.46 x 10 -9 ergs cm~ 2 A -1 (marked with an 
"x" on the figure). The scaling factor of 6.44 x 10 -17 , 
which, in the absence of interstellar extinction, is es- 
sentially the square of the angular radius, implies an 
angular diameter for Vega of 3.31 mas. This is rea- 
sonably consis_te^rt_with_^^ of 3.24±0.07 
mas bv iHanburv Brown et al.l (|1974f ). 3.28±0.01 mas by 
iCiardi et al. (2001), and 3.32 mas bv lAufdenberg et al.1 
(|2006| ). the last of which was based on rapidly rotating 
models of Vega with a large equator-to-pole temperature 
range. As can be seen in Figure [5j the model provides 
an excellent overall representation of the SED, particu- 
larly in the relative levels of the optical and UV fluxes. 
The major discrepancy is that the Balmer continuum in 
the model appears less steep than in the observations. 
This results in an overestimate of the flux just short- 
ward of the Balmer jump and an underestimate near 
~2000 A, both by no more than 5-6%. The observed 
and predicted fluxes are consistent between about 1300 
and 1500 A. Also, although not obvious in Figure [5j the 
model overestimates the flux between the upper Balmer 
lines and produces somewhat narrower Balmer lines than 
seen in the observed SED. A higher surface gravity, near 
log g = 4.0, would resolve the Balmer line width discrep- 
ancy, but not the Balmer continuum problem. 

In Figure [6] we zoom in to provide a closer view of the 
comparison between the model and observed SEDs in 
the UV. The Balmer continuum discrepancy is evident in 
the top panel, which illustrates the wavelengths covered 
by IUE's LWP camera, although the detailed structure 
in the model SED matches the observations beautifully. 
The middle panel, covering the range of the SWP camera 
also shows good detailed agreement between model and 
observations, although the model fluxes tend to run high 
longward of 1600 A. Recall that the A < 2000 A region is 
where our processed low-resolution IUE SWP spectrum 
disagrees with the CALSPEC SED and we cannot rule 
out that some of the differences between model and ob- 
served SEDs are due to calibration issues. Nevertheless, 
the best-fit model reproduces the UV SED level (as nor- 
malized by the flux at 5556 A) very well and we conclude 
that results based on the UV line spectrum and the UV 
SED are in accord. 

The bottom panel in Figure [6] shows the small island 
of flux below Lyman-a. This region was not included 
in the fitting procedure and it is clear that the model 
fails to reproduce the observed flux level as revealed by 
the BEFS data. The strong absorption features observed 
in the model and observed spectra are ground state C I 
lines. The discrepancy cannot be remedied simply by in- 
voking a larger value of T c g , since this would require a 
temperature increase of much more than 100 K (which 
is clearly inconsistent with the rest of the spectrum) and 
would not match the shape of the BEFS energy distri- 
bution. The shape of the emergent flux distribution in 
this region is modulated by the Lyman-a line on the red 
side and the confluence of strong C I lines on the blue, 
culminating in the ground state C I ionization edge at 
1102 A (see, e.g.. lLemke k, Vennl[l996l ). In addition, the 
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region is sensitive to the temperature structure of the 
outer atmosphere and shows weak emission features of 
Fe II and Cr II that cannot be accounted fo r by LTE cal- 
culations (see, e.g.. Ivan Noort et~aT1 119981 ). The NLTE 
calculations of Ivan Noort et al.l are also able to account 
for the non-zero flux in the Lyman-a line core, as seen in 
Figure |3~T1 Even in the absence of Vega's own particular 
peculiarities, it is perhaps not surprising that our LTE 
model predictions breakdown in this spectral region. 

We conclude that our best-fitting atmosphere model 
for Vega — determined exclusively from fitting the ab- 
sorption line spectrum in the UV from 1282-3097 A and 
without reference to the shape or level of the observed 
SED — produces excellent agreement with the emergent 
flux spectrum over the range ~ 1200-7000 A. The chief 
discrepancies noted are the slope of the Balmer contin- 
uum, the overestimate of flux between the higher or- 
der Balmer lines, and the underestimate of the flux in 
the 1100-1200 A spectral island. It is not immediately 
clear if these discrepancies can be attributed to the gen- 
eral physical and computational complexities inherent in 
modeling these spectral regions — and thus would limit 
our ability to model the spectra of any A-type star — or 
if they represent a particular problem affecting Vega due 
to its own physical complexity. 

6. THE COMPOSITION OF VEGA'S ATMOSPHERE 

Given the beautiful agreement between the observed 
and computed spectra for Vega, and the consistency be- 
tween spectral line and SED diagnostics, we proceed with 
some confidence to examine the composition of Vega's 
atmosphere are revealed by our best-fit model. The un- 
certainties listed in Table [2] for the 17 elements with in- 
dividually measured abundances reflect only — as noted 
earlier for the other parameters — the statistical uncer- 
tainties in the measurement technique. Additional un- 
certainty comes from some of the input assumptions in 
the analysis and we have attempted to gauge the im- 
pact of these by running the other fits listed in Table [2j 
In FITS B and C, we changed the assumed abundance 
of the Kurucz atmospheric structure model (and of all 
elements not determined individually) by ±0.5 dex, re- 
spectively and in FITS D and E we changed the assumed 
value of i>sin« by ±3.0 kms -1 , respectively. To further 
investigate the stability of the abundances, we also forced 
changes of ~ ±0.2 dex in the best-fit value of the surface 
gravity \ogg (FITS F and G) and about ±0.5 kms -1 in 
the best-fit value of the microturbulence velocity Vturb 
(FITS H and I). These latter two sets of tests are rather 
heavy-handed in that they deliberately skew the fit away 
from the optimal values, but nevertheless provide some 
insight into its robustness. 

All in all, as a detailed look at Table [2] will reveal, the 
abundance determinations are remarkably stable. Par- 
ticularly note the virtual insensitivity to the choice of 
v sin i and the small variations induced by a factor-of- 
10 change in the composition of the ATLAS9 structure 
model. This latter point is important since the scaled so- 
lar composition of the ATLAS9 models used for our anal- 
ysis does not match the composition determined from 
the UV spectral synthesis. We used the [m/H] = —0.5 
models since they are the closest ATLAS9 models to the 
overall composition inferred from our results, but the de- 
tailed element-by-element composition profiles are differ- 



ent. An obvious follow-up to our study would be to redo 
the fitting procedure using a grid of atmosphere models 
computed with a more Vega-like composition, e.g., using 
the A TLAS 12 program (|Kurucdll99"alCastelh fe Kuruczl 
Il994h . However, the insensitivity exhibited by our anal- 
ysis to large changes in [m/H] suggests that this would 
have li ttle effect on our re s ults, consistent with the find- 
ings of iCastelli fe Kuruczl ([1994). The largest variations 
in the abundances (of some elements) occur, not sur- 
prisingly, when the microturbulence velocity is disturbed 
from its best fit-value. The particular elements most af- 
fected should be those whose most important lines (from 
the point of view of the abundance determination) lie of 
the "flat part" of the curve of growth and are thus most 
sensitive to changes in v tur b- Stronger and weaker lines 
are not affected. Table [2] reveals that the elements S, Fe, 
Ni, and Zn are most affected by ad hoc modifications in 
vturb, which indicates that they are most influential in 
determining its value. Considering all the tests in Table 
[2l it is not unreasonable to consider the abundance deter- 
minations to be generally stable against the assumptions 
in the analysis to a level of ~ ±0.05 — 0.07. 

Other fundamental sources of uncertainty in the abun- 
dance results, not accounted for in Table [2j arise from 
concerns about the available atomic data and the input 
physics used in the models. It is in this regime that 
the benefits of our approach are most obvious. Namely, 
the use of observations in the UV spectral region and 
the adoption of the latest atomic transition data mean 
that, for a given element, the analysis generally utilizes 
spectral features from low-lying atomic levels, includes 
dominant ionization stages and sometimes more than one 
ionization stage, and potentially utilizes many different 
lines arising from a number of different energy levels. 
These attributes reduce the sensitivity of the results to: 
1) NLTE effects (which are generally reduced for low- 
lying levels, see below), 2) com plications due to t he rota- 
tional broadening formulation (|Yoon et"ai1l2008l ). and 3) 
uncertainties in the atomic parameters for specific tran- 
sitions. 

Before examining the composition of Vega's atmo- 
sphere, we first comment on the abundance determina- 
tions for several specific elements: 

• CARBON The abundance of C is particularly sta- 
ble and is based on numerous C I lines and two 
strong C II lines near 1335 A (see Figure 13. 3|) . 
The C I and C II results are compatible, i.e., 
excluding the C II lines from the analysis pro- 
duces nearly no change in the derived C abun- 
dance. When first using the original log gf val- 
ues obtained from VALD, we saw that C I lines 
from the 10192.63 cm -1 level produced large dis- 
crepancies between the model and the observations. 
NLTE effects are not expected to be strong for 
this low-lying level (|Stuerenburg fc H olwcgcr 1990) 
and poorly-determined atomic parameters seemed 
a likely source for the discrepancy, which varies in 
severity among the different lines. Lines from this 
level were excluded from the analysis by setting the 
data weights to zero in their vicinity and empiri- 
cal estimates of their log gf values were made, as 
detailed in Paper II. The fits shown in Figures 12.11 
13.21 and !3.3l were performed with these empirically- 
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determined values although, again, these features 
did not contribute to the determination of the C 
abundance. The ground state C I lines (i.e., from 
the 0, 16.4, and 43.4 cm -1 levels) are well-fit by 
our model, as has been seen in Figures 13.21 and 
13.31 and also in the top panels of Figure [7] where 
we show closeups of the 1560 A and 1657 A com- 
plex es. These panels can be compared with Figure 
9 of IGarcia-Gil et all (|2005| ). Our results suggest 
th at the difficulty in fittin g those lines encountered 
by lGarcia-Gil et aLl (|2005| ) was due to the adoption 
of a too-high C abundance, and was not an NLTE 
effect. 

NITROGEN The N abundance is also very stable, 
as can be seen in Table [5] The results are based 
on a number of strong N I lines from the 2s 2 2p 3 
levels at 19224, 19233, 22839 cm" 1 . The good fit to 
the strongest N I lines, at 1492.625, 1492.820, and 
1494.675 A, can be seen in Figure [7] The resultant 
abundance is also compatible with — although was 
not based on — the resonance triplet near 1200 A, 
as can be seen in Figure 13.11 NLTE effects are not 
expected t o be significant for t he low-lying levels 
used here (|Rentzsch-HolmllT996ft . 

OXYGEN The O abundance is strongly influenced 
by the very strong lines at 1302.168, 1304.858, and 
1306.029 A, arising from levels at 0, 158, and 227 
cm -2 , respectively. These lines can be seen in Fig- 
ure 13.21 The latter two lines are blended with 
strong Si II features at 1304.370 and 1305.592 A. 
Much weaker O I lines at 1355.598 and 1358.512 
A are coincident with C I transitions from the 
10192.63 cm -1 level and are masked from the anal- 
ysis. 

MAGNESIUM The Mg abundance is based on a 
large number of Mg I and Mg II lines. The most 
prominent Mg features are Mg II 2797.929 and 
2799.998 A and the Mg II h and k resonance lines at 
2795.528 and 2802.705 A. The region surrounding 
these features is shown in Figure[7] The broadening 
of the h and k lines was handled by SPECTRUM 
using the Anstee-O'Mara theory of line broadening 
by co llisions with neutral H atom s, with coefficients 
from lBarklem fc O'Maral (fl99l ). 

SILICON The Si abundance is based on numer- 
ous lines of Si I and Si II, including strong Si II 
resonance and fine structure features at 1304.370 
and 1309.276 A (see Figures U and O and at 
1526.707 and 1533.431 A (see the closeup view in 
Figure [7]). The derived abundance produces good 
agreement with the strong resonance features near 
1195 A ( see Figure 13. ip , although this wavelength 
range was not included in the fit. Our original line 
list produced large discrepancies between model 
and observations for a number of Si I features, 
mainly between 1539 and 1563 A, arising from low- 
lying levels. As in the case of C I, these levels 
are not expected to be susceptible to NLTE ef- 
fects and the varying level of the discrepancy from 



feature-to-feature suggests deficiencies in the tran- 
sition strengths. We empirically measured log gf 
values for many of these features (see Paper II) 
and utilized these results in producing the spectra 
seen in Figure 2 — although the lines themselves 
did not contribute to the Si abundance. 

• PHOSPHORUS The element P stands out in our 
analysis as showing the largest under- abundance 
relative to the solar standard (see below). The 
abundance is based on a number of weak lines aris- 
ing from the ground state and low-lying levels of 
P I and P II, all located between 1301 and 1787 
A. P has one of the largest statistical uncertain- 
ties among the abundance determinations — re- 
flective of the weakness of the lines and their fre- 
quent blending with other features — but not large 
enough to explain its apparent depletion relative to 
other elements. 

• ZINC The element Zn stands out in our analy- 
sis as showing the only over-abundance relative to 
the solar standard (see below). The abundance 
is strongly influenced by the two strong Zn II 
resonance lines at 2025.483 and 2062.005 A. The 
stronger of these is shown in closeup in the bot- 
tom row of Figure [7] Both features are located 
in noisy regions of the spectrum and are blended 
with nearby lines, resulting in the largest statis- 
tical uncertainty among the abundance measure- 
ments. There is no reason to question the accu- 
racy of the oscillator strengths for the these tran- 
sitions since interstellar medium studies utilizing 
these lines typically find the Zn abundan ce to be 
comp arable to the solar standard (e.g., Uenkinsl 
[2009h . 

Given the above background, we present in Figure 
[8] the composition results from our analysis of Vega's 
UV spectrum. The top panel shows the abundances 
(filled circles) in the standard form relative to the abun- 
dance of H, i.e., as listed i n Tab le [2] Open circles 
show the IGrevesse fc Sauval (|1998L hereafter GS98) so- 
lar composition for comparison. The abundance val- 
ues are plotted against atomic number and the ele- 
ments are identified along the top of the figure. The 
bottom panel of Figure [5] plots the difference be- 
tween the Vega and the GS98 solar composition in 
the form [X/H] = log (N x /N H ) yesQ - log (N X /N H ) G598 . 
The small x's in the panel show the results from 
other investigations of Vega's atmospheric c o mpos i- 
tion, a nd in cl ude re s ults f ro m iFerrero et al.l 1 1983). 
IGieasI (119861) . IGieasI (11988ft. iStuerenbure fc Holweeer 
(11990ft . HAdelman fc Gulliver! (11990ft. iLemke fc Venn 



( 1996 ?) . iRentzsch-Holm 
iWedemeverl 



.19961) llliiic et al.l 
. iu et all (12001ft. iPrzvbilla et al.l 
(120011). ISaffe fc Levatol (|2004ft . IGarcia-Gil et al.l ffixW . 
and lYoon et all (|2008f) . 

Our results are not sufficiently different from the gen- 
erally accepted pattern of abundances to warrant a de- 
tailed discussion of Vega's composition, but a number of 
significant points can be made: 

1 . Figure [8] shows that the IUE high- resolution data 
and our analysis technique allow access to an im- 
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portant range of elements, including CNO, the light 
metals (e.g, Mg, Si, etc.) and the Fe-group. 

2. The established abundance pattern of Vega (i.e., 
significant depletion of the elements heavier 
than Mg and more s olar-li ke CNO), which lead 
iBaschek fe Slettebakl (|1988l ) to note its resem- 
blance to the A Boo stars, is evident in our results. 

3. The depletions of Mg and heavier elements show 
good element-to-element coherence, with the ex- 
ceptions of P and Zn, as discussed above. The 
mean depletion of these elements (excluding P and 
Zn) is [X/H] = -0.72 ± 0.13 (s.d.), i.e., a factor of 
five subsolar as compared to the GS98 solar scale. 
The scatter is larger than can be accounted for by 
the uncertainties in the analysis. 

4. The slightly super-solar abundance of Zn is not 
inconsistent with the hypothesis that the A Boo 
(and Vega) abundance phenomenon arises from 
the accretion of metal-deficient circumstellar gas, 
with the deficiency arising from t he condensation of 
some elements onto dust grains ()Venn fc Lambert] 
1990). In the interstellar medium, Zn is seen to re- 
sist incorporation into grains and generally mai n- 
tains a solar-level abundance (e.g.. Uenki ns 2009). 

5. The enhancement of N relative to C and O is con- 
sistent with the hypothesis that processed mate- 
rial may have been mixed to the surface, as a re- 
su lt of Vega's rapid rotation. (See the discussion 
in lYoon et all 120081 ). On the other hand, the de- 
pletion characteris tics of N in th e ISM differ from 
those of C and O (| Jenkins! I2009f ) and accretion of 
depleted circumstellar gas might produce such a 
pattern. Both processes may be occurring simulta- 
neously. 

7. FINAL COMMENTS 

We have shown that the emergent UV spectrum of 
Vega, as recorded by the IUE satellite in both high- 
dispersion and low-dispersion, can be fit remarkably well 

- both in detail and in its broad properties — by 
a single-temperature synthetic spectrum based on the 
LTE, plane-parallel ATLAS9 models of Kurucz and the 
LTE spectral synthesis program SPECTRUM of Gray, 
with an updated atomic line list. The chief discrepancies 
between the model and the observations — i.e., those 
which cannot be eliminated by adjustments in the model 
parameters — are an underestimate of the slope of the 
Balmer continuum, producing a model flux excess near 
3500 A and a deficit near 2000 A (both by up to 5-6%) 
and a deficit of model flux in the spectral "island" be- 
tween 1100 and 1200 A. Since few A-type stars have UV- 
through-optical data approaching the extent and quality 
of Vega's, it is not clear whether these deficiencies are 
Vega-related, and result from its high rotational veloc- 
ity and subsequently large gravity-darkening, or are gen- 
eral deficiencies which would plague modeling attempts 
of any A-type star. Studies of additional, more normal, 
stars are needed to address this issue. At the detailed 
level, small discrepancies at discrete wavelengths point 



to incompleteness and to deficiencies in the atomic line 
list used to generate the model. The most notable dis- 
crepancies involve transitions from the 10193.63 cm -1 
level of C I and ground-state and fine-structure levels of 
Si I. In conjunction with a concurrent study of the UV 
spectrum of Sirius A, we have produced empirical esti- 
mates of the log gf values of these lines (see Paper II) . 
Attempting to improve the quality of the UV transition 
data will be an ongoing project. 

If Vega were a normal, equator-on, slow-rotating star, 
then its emergent spectrum and our analysis — which 
is based strictly on the high-dispersion UV line spec- 
trum, but consistent with the UV-through-optical SED 
— would indicate that it has a surface temperature T ff 
~ 9550 K, surface gravity logg ~ 3.7, a general sur- 
face metallicity [m/H] ~ —0.5, and a microturbulence 
velocity Vt ur b — 2.0 kms" 1 . The Balmer lines, which we 
examined only cursorily, seem to require a larger log g to 
explain their widths, as has been found in other studies. 
Since Vega is not a normal star but, rather, a rapid rota- 
tor, these properties must be treated with caution since 
the emergent spectrum will be dependent on viewing an- 
gle. Our results must be regarded as representing some 
sort of average across the observed hemisphere. In gen- 
eral, the simple spectroscopically-derived properties are 
not the values that one would adopt in stellar structure 
calculations to understand Vega's evolutionary state. 

Modeling the complex UV line spectrum of Veg 
revealed by the IUE high-dispersion data, has allowed us 
to determine the specific surface abundances for 17 differ- 
ent chemical elements, including CNO, the light metals, 
and the iron group elements. The resultant abundance 
pattern agrees in general with previous results, although 
there is considerable scatter in the literature. The suc- 
cess of the modeling technique, coupled with an a prior 
appreciation of the advantages to modeling transitions 
from low-lying levels, transitions from multiple ionization 
states of the same element, and, in some cases, hundreds 
or more transitions per element, suggest that these are 
the most reliable measurements to date for the respective 
elements. 

Despite its peculiarities, Vega has turned out to pro- 
vide a powerful test of the extent of our abilities to model 
the atmospheric properties of the early A-type stars, par- 
ticularly the detailed UV line spectrum. The value of the 
measurements from this pilot study will increase as this 
analysis is extended to more objects in the rich high- 
dispersion IUE data archive, including both normal and 
peculiar objects. 

The author acknowledges support from NASA grant 
NNX08AJ62G. Special thanks to Richard Gray of Ap- 
palachian State University for making his SPECTRUM 
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TABLE 1 

High-Resolution UV Spectra of Vega 



Instrument 


Observation 


Program 


Wavelength 


Spectral 




ID 


ID 


Range 


Resolution 


WE/bWP 


SWP32870 


PHCAL 


1150-2000 A 


23 kms 1 


TTT TP /CTS ^TD 

WE/bWP 


b WP32887 


PHCAL 


1150— zUUU A 


16 km s 


TTT TP fQWT'D 

1 t/ii/bWr 




D 1 J HP 


i 1 cn onnn A 
lloU— zUUU A 


HQ 1,-Kvi n — 1 

ZO Kill S 


TTT TP ICWT'D 

1 t/ii/aWr 


OVVP38411 




i 1 cn onnn X 
lloU— zUUU A 


zo km s 


TT TTP IGWTT3 

1 UM//b\vir 


bWP42o21 




i 1 cn onnn X 
lloU— zUUU A 


23 kms -1 


TTT TP ICWT'D 

1 t/ii/aWr 


bWP45z8o 




t 1 cn onnn X 
lloU— zUUU A 


zo km s 


TTT TP IT WTYi 

1 U-Cj/ LW r 


LWPO/490 


QTUD ID 


i nnn oonn X 

iyuu— ozUU A 


16 kms - * 


TTT TP IT WIT} 

1 U-Cj/ -LW r 


T 1X713 1 OKI 7 


DUP A T 


i nnn oonn X 


16 kms -1 


IUE/LWP 


LWP 17648 


PHCAL 


1900-3200 A 


16 kms -1 


IUE/LWP 


LWP21295 


PHCAL 


1900-3200 A 


16 kms -1 


IUE/LWP 


LWP23642 


PHCAL 


1900-3200 A 


16 kms -1 


#ST/GHRS/G160M 


Z3G05208T 


6828 


1186-1122 A 


16 kms -1 


//ST/GHRS/G160M 


Z3G05206T 


6828 


1276-1311 A 


16 kms -1 


#ST/GHRS/G160M 


Z3G05207T 


6828 


1305-1340 A 


16 kms -1 


JJST/GHRS/G200M 


Z2MW0207T 


5673 


1840-1879 A 


13 kms -1 


#ST/GHRS/ECH-B 


Z1210306T 


2461 


2590-2603 A 


3 kms -1 


tfST/GHRS/ECH-B 


Z121030AM 


2461 


2339-2352 A 


3 kms -1 


tfST/GHRS/ECH-B 


Z121030DM 


2461 


2847-2861 A 


3 kms -1 


BEFS a 


BEFS2094 




900-1220 A 


~60 kms -1 



a The Berkeley Extreme and Far-Ult raviolet Spectrome ter (BEFS) was flown aboard the ORPHEUS telescope on the ORPHEUS-SPAS space 
shuttle missions in 1993 and 1996 (see|D et al.|[2003) 



TABLE 2 

Spectral Synthesis Results for Vega 



Property 


FIT A 


FIT B 


FIT C 


FIT D 


FIT E 


FIT F 


FIT G FIT H 


FIT I 


No. of 




















Lines a 


Tef} (K) 


9547 ± 17 


9546 


9546 


9543 


9535 


9448 


9549 


9545 


9541 




log 9 


3.72 ± 0.03 


3.69 


3.83 


3.71 


3.71 


3.60 


4.00 


3.73 


3.73 




[m/H] b 


-0.50 


-1.00 


0.00 


-0.50 


-0.50 


-0.50 


-0.50 


-0.50 


-0.50 




v turb (km s" 1 ) 


2.04 ± 0.02 


2.06 


1.94 


1.80 


2.21 


2.00 


2.00 


1.50 


2.50 




usini (km s — 1 ) 


21.80 


21.80 


21.80 


18.80 


24.80 


21.80 


21.80 


21.80 


21.80 




v rad (SWP) c (km s" 1 ) 


-16.14 ± 2.46 


-16.10 


-16.12 


-16.03 


-16.18 


-16.12 


-16.09 


-16.16 


-16.05 




iw(LWP) d (km s- 1 ) 
C/H e 


1.37 ± 1.88 


1.35 


1.38 


1.45 


1.32 


1.36 


1.36 


1.32 


1.41 




8.18 ± 0.02 


8.15 


8.22 


8.18 


8.18 


8.14 


8.11 


8.19 


8.16 


196 


N/H c 


7.91 ± 0.03 


7.88 


7.95 


7.91 


7.92 


7.88 


7.86 


7.97 


7.83 


46 


O/H 6 


8.48 ± 0.03 


8.45 


8.52 


8.47 


8.49 


8.44 


8.47 


8.48 


8.49 


8 


Mg/H c 


7.10 ± 0.02 


7.07 


7.12 


7.08 


7.09 


7.03 


7.11 


7.10 


7.08 


327 


Al/H c 


5.57 ± 0.02 


5.55 


5.61 


5.55 


5.57 


5.52 


5.62 


5.60 


5.54 


84 


Si/H c 
P/H c 


6.71 ± 0.01 


6.69 


6.74 


6.70 


6.72 


6.68 


6.68 


6.75 


6.67 


388 


4.21 ± 0.05 


4.18 


4.26 


4.19 


4.22 


4.15 


4.18 


4.23 


4.19 


78 


S/H e 


6.61 ± 0.04 


6.56 


6.65 


6.60 


6.62 


6.57 


6.51 


6.71 


6.47 


75 


Ca/H c 


5.73 ± 0.03 


5.71 


5.78 


5.74 


5.75 


5.71 


5.68 


5.83 


5.67 


91 


Ti/H c 
V/H c 


4.21 ± 0.02 


4.18 


4.27 


4.20 


4.23 


4.16 


4.27 


4.27 


4.17 


514 


3.40 ± 0.03 


3.37 


3.45 


3.36 


3.42 


3.34 


3.44 


3.42 


3.38 


462 


Cr/H e 


5.08 ± 0.01 


5.06 


5.12 


5.06 


5.09 


5.04 


5.13 


5.15 


5.01 


2933 


Mn/H c 
Fe/H c 


4.68 ± 0.02 


4.65 


4.72 


4.65 


4.69 


4.64 


4.73 


4.73 


4.63 


1690 


6.78 ± 0.01 


6.76 


6.82 


6.77 


6.78 


6.74 


6.83 


6.91 


6.65 


7416 


Co/H c 


4.14 ± 0.03 


4.11 


4.19 


4.13 


4.14 


4.09 


4.19 


4.21 


4.08 


637 


Ni/H c 
Zn/H c 


5.37 ± 0.02 


5.36 


5.44 


5.37 


5.38 


5.34 


5.44 


5.53 


5.25 


738 


4.81 ± 0.10 


4.80 


4.87 


4.83 


4.82 


4.78 


4.86 


4.93 


4.72 


12 



Note. — The bold-faced model in the second column (FIT A) is the "best-fit' model discussed in the text. The other fits are described in £|4land 
Ef6l The underlined parameters were held fixed at the indicated values in the various analyses; all others were determined by the \ 2 minimization 
procedure. 

a The approximate number of spectral features contributing to the abundance measurement for each element. 

The quantity [X/H] refers to the general metallicity of the Kurucz ATLAS9 models on which the spectral synthesis calculations were based. 
c The quantity iVaeKSWP) refers to the mean value for the velocity offset between the model spectrum and the 38 IUE SWP orders. The uncertainty 
listed is the standard deviation of the sample. 

^ The quantity V ra d (LWP) refers to the mean value for the velocity offset between the model spectrum and the 42 IUE LWP orders. The uncertainty 
listed is the standard deviation of the sample. 

° The individual elemental abundances determined from the spectral synthesis analyses expressed in the form X/H = log(Nx/Nn) + 12.0, where 
Nx/Nh is the number abundance of clement X relative to that of hydrogen. 
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Fig. 1. — Upper Panel: High-dispersion IUE spectrum of Vega, constructed as described in [[2] from the six short-wavelength SWP 
spectra (A < 1980 A) and 5 long- wavelength LWP spectra (A > 1980 A) listed in Table [T] The apparent noisiness of the data is due to 
the thousands of individual absorption lines, although the short wavelength end of the LWP data (1980-2300 A) has a significantly lower 
signal-to-noise ratio than other regions. Small gaps at the long wavelength ends of the SWP and LWP data arise from incomplete overlap 
between adjacent echelle orders. Some medium frequency structure in the data, most noticeable around 2200 A, arises from incomplete 
removal of the echelle "ripple." The spectral resolution of the SWP and LWP data are ~23 kms -1 and ~16 kms -1 , respectively. Lower 
Panel: Best-fitting model of Vega's UV line spectrum, as described in [0 broadened to the resolution of the IUE high-dispersion data. 
Apart from the many individual absorption lines, the most notable features in the spectrum are the 1400 A Lyman-a satellite feature, a 
C I edge at 1446 A due to ionization from the 21648.0 cm -1 level, a Si I edge at 1520 A due to ionization from the ground state, and the 
strong Mg II h and k lines at 2800 A. 
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1284 



1286 



1288 



1290 



1292 



1294 



x 



o 

Z 




1296 



1298 



1300 



1302 



1304 



1306 



x 
3 



- 

o 

Z 



1.2 




1 


1.0 




i 1 A 1 - 

i a LA . A IT\ k.i k 


0.8 




y^n/ At 




0.6 


- a 1X1 1 ft a sn 


1 V 1 /\Vn 
'11 1 




0.4 


t If \// V/v/^\ 11 / I 111 V 

IT \ MJn A 1 ■ \\ 






0.2 


- SWP order 105 \ / 1/ U\ / 




0.0 
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Wavelength (A) 
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Fig. 2.1. — A detailed view of the model fit to Vega's high-resolution IUE spectrum. Shown arc the SWP echelle orders 107, 106, and 
105, covering the wavelength ranges 1282.22-1294.22 A, 1294.24-1306.49 A, and 1306.49-1318.95 A, respectively. The model and the data 
are indicated by the thick and thin curves, respectively. Data points excluded from the analysis due to detector artifacts are identified with 
an "x," as at 1292.5 A. Figures 2.1-2.27, which cover the full wavelength range and all 80 SWP and LWP echelle orders included in this 
analysis, are available in the online version of the Journal. An ascii file which identifies the various absorption features can be obtained 
from the author. 
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Fig. 3.1. — A comparison between the best-fit model for Vega (thick curve) and GHRS/G160M spectrum Z3G05208T (thin curve). The 
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Fig. 5. — A comparison of the observed UV-through-optical SED for Vega (thin curve) and the best-fit model (thick curve). The model 
was scaled by a factor of 6.44 X 10 -17 (corresponding to an angular diameter of 3.31 mas) to match the observations at 5556 A (marked with 
an "x"). The observed SED was taken from the CALSPEC database, with some modifications in the far-UV (see jj5}. For this comparison, 
the model was extended outside the wavelength range over which the fitting procedure operated (1282-3097 A) and was broadened to 
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Fig. Online 2.1. — A detailed view of the model fit to Vega's high-resolution IUE spectrum, for SWP echelle orders 107, 106, and 105, 
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Fig. Online 2.3. — Same as Figure [Online 2TTI but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 
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Fig. Online 2.4. — Same as Figure [Online 2~T\ but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 
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Fig. Online 2.5. — Same as Figure [Online 2~T\ but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 
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Fig. Online 2.7. — Same as Figure [Online 2TTI but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 
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Fig. Online 2.8. — Same as Figure [Online 24*1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 
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Fig. Online 2.9. — Same as Figure [Online 2TTI but for an additional set of high-dispersion IUE echelle orders, identified in the lower left 
of each panel. 



UV Spectrum of Vega 



31 




0.0 

i7ir 

2.0 



1720 



1725 



1730 



1735 



x 

£ 1.5 

N 1.0 



g 0.5 



0.0 

































V fV * \(f v 1 












VP order 79 
























1740 


1745 


1750 


1755 








1760 



1765 



1770 
Wavelength (A) 



1775 



1780 



Fig. Online 2.10. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.11. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.12. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.13. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.14. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 




Fig. Online 2.15. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.16. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.17. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.18. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.19. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.20. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.21. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.22. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.23. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.24. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.25. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.26. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 
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Fig. Online 2.27. — Same as Figure [Online 2TT1 but for an additional set of high-dispersion IUE echelle orders, identified in the lower 
left of each panel. 



